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Abstract 

In this paper we present a novel project-and-lift approach to compute the set of minimal 
generators of the semigroup (Afll™,+) for lattices A C IT 1 . This problem class includes 
the computation of Hilbert bases of cones {z : Az = 0, z £ R™} for integer matrices A. A 
similar approach can be used to compute only the extreme rays of such cones. Finally, some 
combinatorial applications and computational experience are presented. 

1 Introduction 

Lattice points in polyhedral cones arise as interesting objects in many branches of mathematics as, 
for example, combinatorics, integer programming, computational algebra, or topology. Often one 
is interested in a finite subset of lattice points in the cone that generate all other lattice points in 
the cone as non-negative integer linear combinations. 

Definition 1.1 Let C C M™ be a polyhedral cone with rational generators and let A C Z™ be 
a lattice. Then we call a finite set H = {hi, . . . ,h t } C AflC a generating set of the monoid 
(A (~l C, +) if for every z £ A (~l C there are non-negative integral multipliers Ai, . . . , At such that 

z = yij— i Aj/ij. 

Note that for every pointed rational polyhedral cone C C 1" and every lattice ACZ" there exists 
a unique generating set of the monoid (A n C, +) that is minimal with respect to inclusion (|l2|, 
Chapter 7) . For A = Z" this generating set is also called a Hilbert basis Jl3| . It is in general a hard 
problem to compute Hilbert bases. In this paper we will give a novel project-and-lift algorithm 
to compute minimal generating sets for C = R" and arbitrary lattices A C Z". One prominent 
special case, for which many algorithms have been proposed (e. g. ||, ||, [H], |l2|, 0), is the 
case of A = kerz(A) for matrices A e Z dxn , where ker z (A) = {z : Az = 0, z G Z n } denotes the 
integral kernel of A. In these algorithms, the Hilbert basis is extracted from a (usually much) 
bigger superset, in some cases even additional variables have to be introduced. The advantage of 
the algorithm that we propose here is that we do not need additional variables, that we may throw 
away unnecessary vectors in early and intermediate stages of the algorithm, and that we therefore 
arrive at a far smaller superset from which the set of minimal generators of the monoid can be 
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extracted. A similar project-and-lift approach leads to a novel algorithm to compute the extremal 
rays of cones given as {z : Az — 0, z € R™ }, see Section ||. Note that one may apply this approach 
also for the computation of truncated Graver bases which improves the algorithm given in [ || . 
Some implementations of algorithms to compute extreme rays of cones (e. g. cdd porta ) and 
to compute Hilbert bases (e. g. NORMALIZ are freely available over the internet. The latter 
code, however, needs a generating set of the cone {z : Az — 0, z £ R™ } as input. Thus, no direct 
computational comparison between NORMALIZ and the algorithm presented below is possible. 
The main objective of this paper is an algorithmic solution to the following problem. 

Problem 1.2 Given {p\, ... ,p s } C U 1 that generate a lattice A C Z™ over Z. compute the unique 
minimal generating set H(A) of the monoid (A n R™, +) and the set i?(A) of extremal rays of the 
pointed rational cone Ca = {z : z = 531=1 MPh ^li ■ • ■ >^s G R, z G R™ }. 

In Section || we present a selection strategy for the S- vector that has to be considered next in 
the Completion Algorithm |2.3| to compute the minimal generators of the monoid. This strategy 
leads to a tremendous speed-up of our project-and-lift algorithm. In Section || we collect some 
applications of Hilbert bases and extreme ray computations. These include, for example, the 
computation of dual cones. The extreme rays of the dual cone of a cone C, however, are exactly 
the normal vectors of the facets of C. Finally, in Section ^, we report on an implementation, MLP, 
and some computational experience while working on the problem of counting magic arrays p| . 



2 Computation of Minimal Generators 

In order to state the algorithm below we need to introduce some useful notation. 

Definition 2.1 For any m > j let irj : R m — > R J be the projection onto the first j coordinates. 
Moreover, let Kj := {ir™(v) : v G A}, Kf := Kj n (R J + _1 x R+), KJ := Kj n (R^ 1 xL). 

Note that Kj is a sublattice of 1? and that Kj and KJ are intersections of a lattice, ATj, with 
pointed rational cones, that is, (Kj~,+) and (KJ ,+) are monoids. Let Hj~ and HJ denote their 
unique inclusion minimal generating sets. Clearly, H(A) — iJ+ since K+ = A n R™ . 
We may assume that the generators {pi, . . . ,p s } C U 1 of A have the following structure, which is 
clearly achievable by suitable elementary integral row operations and permutations of columns: 
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with pi } i > 0, i = 1, . . . , s. 

In what follows, we will start with — {(pi,i)}, which is easily computed, and compute , 
, . . . , inductively. For each step — > H^ +l we will employ a completion procedure ||. 
The input sets for j > s and j < s differ slightly due to the fact that for j < s the map 7r™ is not 
injective. In particular, there are non-zero vectors dgA with ^(v) = if j < s. 
As C Kj, we know that for every h £ there is at least one vector v G A with h = 7r™(i>). 
Thus, for every h G Hj there is some Zi'gZ such that (h, h') G Kj+\. If there are several choices 
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for h' , which can only happen if j < s, we will choose the smallest non- negative integer hi such 
that (h, h 1 ) G Kj+i to form the input set F below. For j < s we have to add the vectors 7r™ +1 (pj + i) 
and — 7r™ +1 (pj + i) to the input set F. 

Lemma 2.2 The set 

U {(h,h') : {h,h') G K j+1 } if] > s, and 

heH+ 

U {(h,h>):(h,h')eK j+1 }u{n? +1 (p j+1 ),-ir? +1 (p j+1 )} if j < s, 
heHf 

generates K^ +1 U KJ +1 over Z + , that is, every z G Kf +1 U Kj +1 can be written as a non-negative 
integer linear combination of elements of F . 

Proof. Let 7r" +1 (u) G Kj+i U where v e A. Since Hf is a generating set for the monoid 

(Kj~,+), we can write 7r™(w) = 0!j7r™(t>i) for some 7r™(«i) G and a, > 0. Therefore, 
7r" +1 (i>) — X) a i 7r j'+i ( v i) nas zeros m the first j components and must be the zero vector if j > s, 
or an integral multiple of 7r" +1 (pj+i) if j < s. This concludes the proof. □ 

In what follows, dashed variables, like h', at the end of a vector, like (h,h'), always represent 
integer numbers. 

Algorithm 2.3 (Algorithm to Compute H^ +1 Li H~ +1 ) 

U {(h,h'):(h,ti) &K j+1 } if J >s, and 

heHj 

U {(h,h'):(h,h')GK j+1 }u{nJ +1 ( Pj+1 ),-^ +1 (p j+1 )} ifj<a 
heHf 

Output: a set G which contains Hj~ +1 U H~ +1 
G:=F 

C:= U S-vectors(/, ff ) 

f,geG 
while C ^% do 

s := an element in C 

C:=C\{s} 

f := normalForm(s, G) 

iff^O then 

C := CU U S-vcctors(/, g) 

G:=GU{/} 

return G. 

It remains to define the set S-vectors(/, g) and the function normalForm(s, G). 

In case that one is interested only in those elements of H(A) whose components lie below certain 
bounds u G (Z + U {oo})", that is -ff (A) fl {z : z < u, z G Z™ }, these upper bounds u can easily be 
used to speed-up the computation by considering only the following set of S-vectors: 



Input: F 



S-vectors((w, v'), (w, w')) := 



{(v + w, v' + w')} if v'w' < and v + w < u, 
otherwise. 
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This choice will be justified below after the proof of Lemma |2.5| . 

Behind the function normalForm(s, G) there is the following algorithm, wherein u v if and 

only if u« < uW, i = 1, ... ,j, {u^+^l < \v^ +1 ^\, and > 0. Note that in case of 

u Ej+i u i u and v lie in the same orthant of R J+1 . 

Algorithm 2.4 (Normal Form Algorithm) 
Input: a vector s, a set G of vectors 
Output: a normal form of s with respect to G 

while there is some g £ G such that g Ej+i s do 

a := min{Ls«/ 5 «J : i = 1, .. . ,j+ 1, <?« ^ 0} 

s := s — ctg 
return s 



Note that the first j components of s always remain non-negative during the normal form algorithm. 
Moreover, the (j + l)st component of s either remains non-negative or non-positive. Since ||s — 
ag\\i < || s ||i, the normal form algorithm always terminates. 



Lemma 2.5 Algorithm 2.c terminates and returns a set containing Hj + i U Hj +1 . 



Proof. First, let us show termination of the algorithm. For this define for any r6l the symbols 
r+ := max{r, 0} and r~ := max{— r, 0}. Next, consider the sequence G\F = {(gi,g[), . . . } as it is 



generated by Algorithm 2.3. This sequence fulfills (gi^g'A (gk,g'k) whenever i < k. Therefore, 



(9ii {g'i) + : (9i) ) ^ (9k, (g'k) + > (9k) ) whenever i < k. Applying the Gordan-Dickson Lemma (see 
for example S) to the sequence {(gi, {g[) + , (g[) ), ■ ■ ■ } Q ^+ 2 , we conclude that this sequence 
must be finite and thus Algorithm |2.3| terminates. It remains to show correctness of the algorithm. 



By G denote the set that is returned by Algorithm 2.3 and let (z, z') G Hf +1 L)Hj~ +1 . By Lemma |2~ 
and F C G, we can write (z, z') as a finite positive integer linear combination (z, z') = ^ on(vi, v'A 
for some on G Z>o and vectors (vi,v'A G G with < Vi < z for all i. From the set of all such 
positive integer linear combinations ^a^Vi^v'A choose one such that ^aj|i^| * s minimal. Note 
that, by the triangle inequality, X^iKI — W\ with equality if and only if v[ \—i z' for all i, that 
is, if and only if z' and all v[ are either both non-negative or both non-positive. 
If X) a *KI — then we get Qj+i ( z , z ') for all i- Since (z,z') G u Hj+n the 

representation (z, z') = ^ c<i(vi, v[), on G Z>o must be trivial, that is, (z, z') = 1 • (z, z') 6 G, and 
consequently, (z, z') G G and nothing is left to prove. 

Hence, we will assume on the contrary that > l z 'l holds. Therefore, there must exist 

(■Utijfjj), (vi 2 ,v[ 2 ) such that v[ v[ < 0. The sum (iv,^,) + (i>i 2 , ^ ) was reduced to during the 



run of Algorithm 2.2 which gives a integer linear combination (v^ , v[ ) + (i>i 2 , ) = ^ /3fc (tuft , u^) 
for some positive integers and some (wfcjW^,) G G. 

Moreover, (3k(wk,vj' k ) Qj+i (v^jv'j) + (vi 2 ,v' i2 ) for all fc, implying that < Wk < z and that 

£&KI = lE&«4l = K +<\ < KI + KJ 

holds. But then (z,z') = £ /3 fc (w fc , wj.) + (a u - 1)^ , v'^) + (a ia - l)0i 3 , v' ia ) + Y^ijti lt i a a i( v i> v 
contradicts the minimality of We conclude = \ z '\ anc l the claim follows. □ 

Note that if there are upper bounds u on z, the relations + u,- 2 < z < u always holds in the 
above proof. This justifies the above mentioned restriction of necessary S-vectors that need to be 
considered if upper bounds u on the variables are given. 
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3 Improving the Project-and-Lift Algorithm 



In this section we improve Algorithm |2.3| such that we need not reduce the S-vectors from C with 
respect to G. Instead, it suffices to check each S-vector only for reducibility with respect to G. 
Moreover, the set C of S-vectors that still have to be considered for reducibility need not be stored 
in memory explicitly. Both properties speed-up the computation, save computer memory and thus, 



allow us to solve bigger problems as by application of Algorithm 2.3 alone 



The major idea of the improved algorithm is to consider the elements (v,v') in C by increasing 
norm ||u|ji = Y^i=i v ^' '■ To this end, let Gi :— {(v,v r ) G UffZ_i : |H|i = i}. Suppose we 

have already computed Go,Qi, ■ ■ ■ ,Gk, for some k > 1. Then the following observation allows us 
to avoid all the reduction steps in order to compute Gk+i- 

Lemma 3.1 Every vector (z,z') G Kj+i U ^j+i "with \\z\\\ < k can be written as a finite positive 
integer linear combination J2 a i(di> 9i) w ^ a i £ ^>o, ij9iid'i) £ G<k '■= Go U Gi U . . . U Gk, and 
(9i,9i) Ej+i ( z > z ') f° r al1 i- 

Proof. Let (z,z r ) 6 Kf +l U Kj +1 with ||^||i < k. Clearly, by definition of H+ +1 U Hj~ +V (z,z') 
can be written as a finite positive integer linear combination a i(9i> 9i) with £ ^>o, (diidi) £ 
H+ +1 yjHj +1 , and [ 9i ,g'i) Ej+i i z i z ') for all i. But (gi,g'i) Ej+i ( z > z ') implies \\gi\\i < \\z\\i < k, 
which, together with (<?i,<?-) £ Hj+i U HJ~ +1 , gives (</i,<^) £ £<fc and the proof is complete. □ 

Note that it is crucial for the proof of the following corollary that, in case of j < s, we always 
choose the smallest non-negative integer h! such that (h, h') G Kj + i when we form the input set F 



to Algorithm 2.2. Moreover, note that Go = if j > s and that Go = {nf+iiPj+l), ~ 7r j'+i{Pj+i)} if 
j < s. 

Corollary 3.2 Let (z, z') G Kf+i ^ ^j+i> 1 1 1 1 x = k + 1, and \z'\ < pj + ij + i. Then there are 
precisely two possible situations: 

1. There exists a vector (v,v') G G<k with (v,v') C (z,z'). Then we already know in advance 
that normalForm((z, z ), G<k) = 0. 



2. There does not exist a vector (v,v') G G<k with (v,v') C (z,z'). Then (z,z') G G, 



fc+i- 



Proof. First, suppose that there exists (v, v') G G<k with (v, v') C (z,z'). Clearly, if j < s, 
(v,v') Ej'+i (z,z') and |z'| < Pj+ij+i imply (i7,t/) g^£ = Wj+iiPj+i), -^j+iiPj+i)}- Tnus we 



have d / and consequently ||z — v\\i < fc. We conclude by Lemma 3.1 that there are finitely 
many (not necessarily different) elements (gi,g[) G G<k such that (z — v,z' — v') = XXSiifi) and 
{gi-, g'i) Ej+i {z — v,z' — v') Ej+i (z, z') for all i. Choosing (v, v') together with these (gi, g[) in the 
algorithm normalForm we obtain normaiForm((z, z'), G<k) = 0. 

On the other hand, if there does not exist a vector (v,v r ) G G<k with (v,v') Ej+i (z,z'), then 
(z, z') G t/fe+i, that is, the vector (z, z') cannot be written as a sum (v\, v[) + (v2, v' 2 ) with v^) G 
(Kj +1 U -f^J+i) \ {0} and («i,?4) Ej+i (-^z'), i = 1,2. To see this, assume on the contrary that 
such vectors (vi, w() and (112, v' 2 ) do exist. Note that both vectors lie in the same orthant of W +1 
as (z,z'). From [vi,v'^j Ej+i (*z,z'), i = 1,2, and z' < p 7+l j+1 we conclude that Uj ^ 0, i = 1,2 



Therefore, \\Vi\\i < fc, i = 1,2 which implies by Lemma 3.1 that (vi,v[) and (i>2, V2) can both be 
written as finite positive integer linear combinations of elements from G<k that all lie in the same 
orthant as (vi,v[) and (i>2, v' 2 ). Put together, both combinations allow a finite positive integer linear 
representation of (z, z') by elements from G<k that all lie in the same orthant as (z, z'). But all the 
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summands (v,v r ) € G<k in this representation of (z,z') fulfill (v,v r ) Ej+i ( z , z ') in contradiction 
to our initial assumption that no such non-zero vector (v, v') exists. Thus, (z,z') G Gk+i- D 



In view of Corollary 3.2, considering the elements (v,v') in C by increasing norm ||i>||i leads to an 
algorithm to compute H^ +1 U Hj +1 that has the following advantages: 

• The set C of S-vectors need not be stored, as for k = 2,3,..., the elements in Gk+i are 
generated and immediately checked for reducibility within a simple loop. No superfluous 
vectors are stored. 

• Thus, precisely Hj + i U HJ +l is computed. 

• There is no reduction of vectors, only much cheaper reducibility tests are done. 

• Upper bounds on variables can be easily be used directly to truncate the output set and to 
speed-up the computation. 

It remains to give a criterion of when Hj +1 U HJ +1 = G<k that allows us to stop the inductive 
construction of the Gk- 



Lemma 3.3 Let k 6 Z>o satisfy Gk+i = • • • = G^k = Then H^ +1 U H, +1 = G 



t<k- 



Proof. Note that the input set F defined in Section ^ satisfies F C G<k- Therefore, by Lemma 



2.2, also G<k generates Kj +l U K^ +1 over Z+. Thus, by Lemma 2.5, Algorithm 2.3 would return 
a set contain Hj +1 U HJ +l if it were started with G<k as input set. However, by Lemma 3.1 and 
Corollary 3.2, all S-vectors reduce to with respect to G<k and Algorithm 2.3 returns exactly G<k- 



As all elements in G<k are minimal, we must have G<k = U ^j+i- '— ' 

For j > s, one question remains open from a practical point of view: As the remaining variables 
(with index bigger than j) can be freely permuted and then re-permuted at the end of the com- 
putation, what is a good selection strategy to choose the (j + l) st variable to be appended next? 
First computational experiments show that huge differences in sizes of intermediate sets H~!j~ and 
in running times do indeed occur. 



4 Computation of Extremal Rays 

Given a cone C = L n K" , where L is the linear space generated by {pi , . . . , p s } C Z™ . We 
may again assume that these generators have the same structure as in Section^]. Analogously to 
the computation presented in Section ||, let Kj := {7r™(v) : v £ L}, := Kj n x 

Kj := Kj D (M^ 1 xL). Note that Kf and Kj are pointed rational cones. Let and Rj 
denote their sets of extreme rays. Again we will start with = {(pi,i)}, and compute R% 



R£ inductively. For each step R^ — > Rj+i we will again employ Algorithm 2.3. However, 
we have to use different specifications for the input set, the set of S-vectors, and the algorithm 
normalForm. Note that a similar speed-up as presented in Section [| is possible for this algorithm, 
by choosing elements (v,v') from C by increasing value of |supp(w)|. We will, however, skip the 
details here. 

The following Lemma can be proved in a similar way as Lemma |2.2|. 



G 



Lemma 4.1 The set 



F 



U {(h,h'):(h,h')£K j+1 } ifj>s, ana 

h£R+ 

U {(h, h') : (h, h') € K j+1 } U {^ +1 (Pj+i), -*? +1 (p i+ i)} if 3 < s, 



heRJ 



generates K~!f +1 U K- +l over K + , that is, every z G K^ +1 U K- +l can be written as a non-negative 
linear combination of elements of F. 

Let F as defined in Lemma |4.1| be the input set to the completion algorithm and claim the output 
to be a set G which contains a positive scalar multiple for each element in R^ +1 UR~ +1 . It remains 
to define the set S-vectors(/, g) and the function normalForm(s, G). 

S-vectors((,y), {w,w')) := ( «" " ^ W >^ ~ ^»'» *™ < °> 
" ' v " \ otherwise. 

Note that all S-vectors have as their last component. Behind the function normalForm(s, G) 
there is the following algorithm. 

Algorithm 4.2 (Normal Form Algorithm for Extreme Ray Computation) 
Input: a vector s, a set G of vectors 
Output: a normal form of s with respect to G 

while there is some g £ G such that supp(g) C supp(s) do 

a := min{ S W/ ff W :i = l,... ,j,g (i) ^ 0} 

s := s — ag 
return s 

Note that the first j components of s always remain non-negative during the normal form algorithm. 
Moreover, supp(s — ag) C supp(s) and therefore, the normal form algorithm must always terminate. 

Lemma 4.3 Algorithm \2.3[ with the new specifications for the input set, S-vectors and normalForm 
(as presented in this section) terminates and returns a set containing Rj + i U RJ+i- 

Proof. Termination of the algorithm is clear, since only a finite set C of vectors is checked. Thus, 
it remains to show correctness of the algorithm. 

By G denote the set that is returned by the completion procedure and let (z, z') £ R^ +1 URJ +1 . By 
Lemma L2 and FCG, we can write (z, z') as a positive linear combination (z, z') — ^ ai(vi, v[) 
for some on G M>o and vectors {vi^v'A G G with supp(i>i) C supp(z) for all i. From the set of all 
such positive linear combinations ^2 ai(vi,v'A choose one such that * s minimal. In order 

to find such linear combination we have to solve the linear program 



Q>0 Z ' Z J 



for given (vi,v'f) G G and (z,z'). This linear program is bounded from below by 0. Thus the 
minimal value is indeed attained by some choice of the on. 

Note that — \ z '\ with equality if and only if z' and all v[ are either non- negative or 

non-positive. 
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If X) a «KI = then all {vi,v[) lie in the same orthant as {z,z') and we get snpp((vi , v'A) C 
supp((z, z')) for all i. Together with (z,z r ) = Y^ a i( v i) v l)> a i e ^>o> an d {z, z') G Rj+i U RJ+i 
this representation must be trivial, that is, (z, z') — a\ ■ (v±, v[) G G, a.\ > 0, and consequently, G 
contains some positive scalar multiple of (z, z') as claimed. Hence, we will assume on the contrary 
that X^ a *KI > \ z '\ holds. Thus, there must exist ( v i 2 : v i 2 ) sucn that v' i2 < 0. 

The sum (v n ^d) - (u ^ /« [ 2 ) (v i2 , u - 2 ) = (u^ - {v' ix /v' i2 )v i2 , 0) was reduced to during the run 
of Algorithm 2.3 which gives a linear combination (v^ — (t^ /v' i2 )vi 2 , 0) = flk(u>k, w' k ) for some 
positive scalars (3k and some (wk, w' k ) G G. Moreover, supp((wfe, w' k )) C supp((Uj 1 — (u^/t;^)!!^, 0)) 
for all fc, implying that 

= ^/3feK| = |^^ feU; fcl = K - WiJv'iJv'iJ < \v' n \ - (v'ijv'ijlv'ij 

holds. But if we now choose a — mia.{oii 1 ,ai 2 /(—v' i /v > i )} = minjo;^, — ^i 2 v i 2 / v L} > 0, the 
representation 

(z,z') =^2af3 k {w k ,w k ) + (a h - a) (v h , v' tl ) + (a h + a(u l ' i /^ 2 ))(v i3 ,^ 2 ) + ^ a^i,^) 
contradicts the minimality of ^ a »K'l- We conclude X^ibil = \ z '\ an d the claim follows. □ 



5 Applications 

5.1 Decomposition of Vectors 

Problem. Given A G Z dxn and u G Z™ , find a representation u = 53;=i with hi G 
ff(kerz(A)), i = 1, . . . , i, and positive integral multipliers Ai, . . . , At. 

Solution. Compute a truncated Hilbert basis H(keiz(A)) with upper bounds u on the variables. 
Once -ff(kerz(^4)) has been computed, u can be easily decomposed as a non-negative integer linear 
combination of elements of H(kerz(A)) as follows: Find v G H(ker%(A)) such that v < u and 
replace u by u — v. Note that u — v can also be written as a non-negative integer linear combination 
of elements of H(kerz(A)) as u — v < u. 

This solution approach can be sped-up once an element of some is found that extends to a 
vector v G kerz(A) n K™ with v < u. Note that this vector v must belong to as already its 
projection onto the first j components is not decomposable. In case such a vector v has been found, 
we may replace u by u — v and use the stronger upper bounds u — v instead of u. Note that the 
computation of the truncated Hilbert basis need not be started again. It suffices to throw away 
unnecessary vectors from (that do not lie below the new upper bounds u — v) and to continue 
with the computation of Hj+\- If H£ is reached, we can completely decompose the current u or 
prove that it is not further decomposable and therefore u G H (kerz(A)) C H (keiz(A)) . 

5.2 Counting Magic Arrays and t-designs 

Problem. Given a d dimensional n x n x • ■ • x n array with n d non-negative integer entries, 
wherein the sum of the entries in any coordinate (and in the main diagonal) direction equals a 
magic constant s. Find a formula for the number of magic arrays with given magic constant s. 
Start of Solution. The method presented in |l]] computes in a first step all non-decomposable 
magic arrays. This step is equivalent to the computation of the Hilbert basis of the cone of all 
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magic arrays with non-negative real entries. It turned out that the computation of Hilbcrt bases 
is currently not the computational bottleneck in the method described in 
The computation of Hilbert bases is also employed in counting i-designs 

5.3 Integer Programming 
Problem. Given A e Z dxn , b G Z d , c G R T \ consider the inte ger program 

min{c T z : Az = b, z G {0, 1}"}. 

z 

Given a 0-1 solution zq to Az = b, find a cheaper solution or prove optimality of zq. 
Solution. Transform the problem via z^ — > 1 — {z')^ whenever z^ — 1, z^ — > otherwise, 
into a 0-1 problem with problem matrix A! ', cost vector c', and current feasible solution Zq = 0. 
The above problem reduces to the computation of an element v < 1 in the truncated Hilbert basis 
of ker^(A') n MIL (with upper bounds 1) that has a strictly negative objective value (c') J v. 

5.4 Dual cones 

Problem. Given a cone C generated by integral vectors pi, ■ ■ ■ ,p s in R n - Compute the generators 
R(C D ) and the minimal Hilbert basis H(C D ) of the dual cone 

C D := {v G R™ : z T i; > 0,Vz 6 C} = {d £ M n : pju > 0, i = 1, . . . , s}. 

Note the R(C D ) gives exactly the set of normal vectors of the facets of C. 

Solution. Let I s denote the sx s identity matrix and let P be the sxn matrix that has p\, . . . ,pj 
as rows. Moreover, let L := {(v,u) e K' i+S : Pv — I s u = 0} be generated over R by the vectors 
(ei, Pei), . . . , (e n ,Pe n ). Note that these vectors also generate kerz(P\ — I s ) over Z. Therefore, we 
have C D = n (R" x R+)) and thus 

R{C D ) = {»£ R" :u = Pve R{P r )} and H{C D ) = {v e Z" : u = Pv G H(P T )}, 

where i?(P T ) and H(P 1 ) denote the set of extreme rays and the Hilbert basis of the cone which is 
the intersection of M. s + with the subspace of R s spanned over R by the columns of P (rows of P 1 ). 

5.5 Hilbert Bases of Cones given by their Generators 

Problem. Given a cone C generated by integral vectors pi, ■ ■ ■ ,p s in R", compute the minimal 
Hilbert basis H (C) of C. 

Solution. Since C = (C D ) D , this problem can be solved by first computing the generators R(C D ) 
of the dual cone and by then computing the Hilbert basis of the dual of the dual cone C D , which 
is the Hilbert basis of C. 




6 Computational Experience 

The algorithm presented in Sections |^ and ^ were implemented into a computer program that 
we called MLP. It is written in C and downloadable from tittp://www. testsets.de. Besides 
(truncated) Graver test sets ||, MLP allows the computation of a (truncated) set of minimal 
generators of the monoid A n R™ for any lattice A C Z", in particular for A = kerz(A). 
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In MLP was used to solve some non-trivial instances. For example, it took about 10 days 
on a 1GHz PC with 4GB RAM running linux to compute the set of all minimal 6x6 magic 
squares, yielding a Hilbert basis containing 522, 347 elements. These examples can be found on 



the benchmark section of tittp: //www. testsets.de 
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